Giant Quantum Freezing of Tunnel Junctions mediated by Environments 
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We investigate the quantum heat exchange between a nanoj unction and a many-body or electro- 
magnetic environment far from equilibrium. It is shown that the two-temperature energy emission- 
absorption mechanism gives rise to a giant heat flow between the junction and the environment. We 
obtain analytical results for the heat flow in an idealized high impedance environment and perform 
numerical calculations for the general case of interacting electrons and discuss the giant freezing and 
heating effects in the junction under typical experimental conditions. 



Quantum dynamics of tunnel nanoj unctions is gov- 
erned by underlying relaxation mechanisms and nonequi- 
librium effects since even small currents drive a nanoj unc- 
tion well out of equilibrium [H, [I| . At low temperatures 
the direct energy transfer to the phonon bath becomes 
inefficient and relaxation is dominated by the energy ex- 
change between the tunnelling electrons and electromag- 
netic environment and/or to many-body excitations in 
the electrodes (hereafter we will refer to both mechanisms 
as to relaxation via the environment) [3|. In this Letter 
we calculate the heat flow between the tunnelling elec- 
trons and the environment, which control the junction 
dynamics, using a non-pert urbative technique based on 
quantum kinetic equations taking into account far from 
equilibrium effects. We show, in particular, that a regime 
exists in which the interaction with the environment gives 
rise to a giant cooling of the nanoj unction. 

The energy exchange between tunnelling electrons and 
the environment is determined by the emission of envi- 
ronment modes with temperature equal to that of elec- 
trons, T e , and the absorption of environment excita- 
tions carrying the temperature of the thermal bath, T env . 
Moreover, not only temperatures, but also the distribu- 
tions of emitted and absorbed environment modes may 
appear essentially different in the far from equilibrium 
regime T e ±2 V > T env ±2 Ti ea d s , where V is the voltage 
across the junction. Our main finding is that this two- 
temperature emission-absorption mechanism gives rise to 
a giant heat flow between the junction and the environ- 
ment (see inset in Fig. [I]). 

In the case of a resistive environment [electromag- 
netic fluctuations in a cold (hot) resistor shunting 
the tunnel junction] Q, 0, 0, the heat flow is Q — 
T g RtC \tci(RtC /r e )kB(T e — T env ) / 'RtC with compara- 
ble temperatures \T e — T env | <C T e + T env . Here Rt 
and C are the ohmic resistance [6| and capacitance of 
the junction, respectively, and r e is the electron energy 
relaxation time, r e <C 1/T e <C RtC. This result well 
exceeds the flow in the quasi-equilibrium approximation 
Q ±2 kB(T e — T env ) / RtC [4j], where the emitted and 
absorbed modes have the same temperature. The large 



factor T € RtC \ti(RtC / r e ) ^> 1 by which the two results 
differ, reflects the elevated effective number of environ- 
ment excitations emitted by charges tunnelling through 
the nanoj unction in the out of equilibrium regime, see 
Fig.ffl 




FIG. 1. (color online) Illustration of the non-equilibrium heat- 
ing effects in a nanoj unction. The electrons traversing the 
junction absorb external photons (incident wavy lines) and 
emit them leading to heating of the contact. The plots show 
the giant heating effect as a function of the difference of elec- 
tron and environment temperatures in the non-equilibrium 
situation compared to the quasi-equilibrium approximation 
at zero bias voltage. The full non-equilibrium analysis gives 
an at least one order of magnitude more pronounced heating 
effect than for the latter case: max(Q/Qo) > 10. 

Model. — The rate of the heat flow between the tunnel 
junction and the environment is given by [7j: 



e {n £ P{e) - [1 + n e ]P(-e)} p{e)de, (1) 



where P(±s) is the probability density for the tunneling 
charge-carrier to lose [gain] the energy e to [from] the 
environment. The distribution function n £ in Eq. (pQ) can 
be interpreted as the distribution function of electron- 
hole pairs that appear at the junction interface just after 
the tunneling process: the hole in the source lead and 
the electron in the drain, Fig. [2ji). If the distribution 
functions at the electrodes are Fermi functions with equal 
temperatures T e , then n £ = {(e — V)Nb(s — V,T e ) + 
(e + V)N B (e + V,T e )}/2s, with N B (e ± V,T e ) being the 
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more. The spectral function p(oj) is the probability of 
the electron-environment interaction and characterizes 
the particular system under consideration. 

To estimate the magnitude of the heat flow Q we first 
expand the distribution function P(e) in Eq. (pQ), in the 
first order in p(e): 



FIG. 2. (color online) a) Illustration of electron-hole pair gen- 
eration in the tunnel junction, resulting in the distribution 
function n £ [Eq. ([T])] of these pairs (environment), b) Com- 
parison of the distribution functions for T = and T = V/2. 



equilibrium Bose distribution function. For this case the 
effective temperature of the nano junction is defined as 

T eS (V, T e ) = Urn n e (V, T e ) = \ coth (2) 

At low applied voltages, V <C T e , the effective tempera- 
ture of the junction T e ^ coincides with the temperature 
of the leads, T e ^ w T e . In the opposite case of high volt- 
ages, V ^> T e , we obtain T e ^ « V/2. The function p(e) 
in Eq. (pQ) is the weight function for a junction between 
two normal metals, Fig. [TJ and can be calculated for any 
choice of the electron distribution function in the leads, 
resulting in p(e) = Ae/Rt- 

Heat flow. — To calculate Q one has to specify 
the probability density, which can be written in a 
form P(e) = dt exp[J(t) + iet], where the func- 
tion exp[J(t)] reflects the fact that tunneling electrons 
acquire random phases due to interaction with the 
Bosonic environment, represented by a set of oscillators 
with non-equilibrium distribution of modes, N^. The 
quasi-equilibrium situation where the distribution func- 
tions of the environment modes are Bose distributions 
parametrised by equilibrium temperatures was discussed 
in Ref. Q. In general far from equilibrium, J(t) is jjjj]: 



CO 

J(t) =2 [ —p(lj) \N^e iu3t + (1 + ^ out ))e 
J ^ L 



(3) 

The mode distribution, N^, is defined by a kinetic equa- 
tion with scattering integral describing the energy ex- 
change between environment modes and tunnelling elec- 
trons. The terms proportional to and l+A^S° ut ^ cor- 
respond to the absorbed and emitted environment exci- 
tations, respectively. The combination = l + N^ ^ + 
is the kernel of the time-independent contribution 
to J(t) describing the elastic interaction of the tunnelling 
electrons with the environment modes. In equilibrium 
N u reduces to the Bose- function and the functional P(oj) 
recovers the result of Ref. [2| . In Eq. (j3]) , the energy re- 
laxation time r e determines the low energy cut-off, since 
the electrons start to equilibrate on larger time scales, 
i.e. the non-equilibrium description does not hold any 
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de 
2k 



ep(e) {n £ (l + ^ out )) - (1 + n £ )N^} 

(4) 



The expression in Eq(j4]) becomes zero if 
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. If the distribution functions are not equal to each 



other, we can expand with respect to their differ- 
ence. We consider the case where the voltage bias at the 
nano junction is zero but the temperatures of electrons at 
the leads and those that comprise the environment are 
slightly different, T e = T + ST/2 and T env = T — ST/2. 
Thus, n £ = n £ (T + ST/2), 7Vj in) = n £ (T - ST/2), 
7V £ (out) = n £ (T + ST/2), where n £ is the Bose distribu- 
tion function. Using Eq. (j4]) in the first order in small 
parameter ST /T <C 1 we find 



ST- 



de 
2^ 



sp(s)n' £ (T)(l + 6n £ (T)), (5) 



where n' e (T) = dn £ (T) / de. The index is for the quasi- 
equilibrium situation when the temperatures of emitted 
and absorbed environment excitations are equal and 1 for 
the non-equilibrium case (the index 1 is skipped through- 
out this Letter). Since n £ (T) in Eq. (j5j) is always positive, 
the following inequality is valid IQo^l < IQ^I? where 
and refer to the heat flux in quasi-equilibrium 
and in non-equilibrium cases, respectively. The interac- 
tion function p(e) in Eq. (0 quickly decays at frequen- 
cies larger than some characteristic frequency cj max - For 
temperatures T > cj max (quantum regime) we can ap- 
proximate n £ (T) « T/e ^> 1 and find 
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oo Tp(uj)duj 
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IQO^I /^iPHdw 



ln(o; max re) > 1. (6) 



Remarkably, in higher orders with respect to p(e) the 
non-equilibrium heat flow Q differs from the equilibrium 
flow Qo by the same factor (see supplementary material). 
This result holds even for a finite electric current flow- 
ing through the junction. However, in this case we need 
to replace the temperature T by the effective tempera- 
ture T e ff, see Eq. (j2j), of the tunnelling electrons. Thus, 
the heat flow between the junction and the environment 
appears much larger than what the quasi-equilibrium es- 
timates predict. 

Ohmic approximation. — We now turn to the sim- 
plest case, an environment with a very high impedance 
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environment interaction 



FIG. 3. (color online) Typical heat exchange Q in Eq. (pQ) 
of the Ohmic environment with the tunnel junction between 
two normal leads. Q(T eff ,T, V) vs T/T and voltages eV/T 

(To = 30UW.)' We USed ^maxKin = 100, T env /T = 1, 

p(0) = 10. Q is measured in units of 10 3 w^ax / (e 2 Rt)- 



as compared to the quantum resistance, Rq. In this 
limit tunnelling electrons easily excite the environment 
modes. The spectral density p(uo) of these modes is 
sharply peaked at the zero frequency, u = 0. For the 
correlation function J(t) in Eq. (j3j) the concentration of 
the environment modes at low frequencies implies that 
the expansion of J(t) over t up to the second order yields 
J(t) « —iat — (b/2)t 2 , where the coefficients a and b 
are defined as a = J^i(l + iV"i, out ^ — N^)p{uj)duj and 
6 = J^i ujp(uj)B UJ duj. Using this expansion for J(t) we 
obtain the following result for the density function P(uS) 



P(s) 



-{e-a) 2 /2b\. 



(7) 



Here the expansion parameter a can be estimated as 

follows a — do ( 1 + (^e-Tenv) ln(u; max T e ) \ w l iere aQ _ 

V TT^max y 

2 J pdcj & 2p(0)co> max ~ 2T? C with E c being the charging 
energy of the tunnel junction, T e is the electron temper- 
ature in the junction, T env is the temperature of environ- 
mental modes, co> max ~ 1/(RtC). Similar for coefficient 
b in Eq. (O we obtain 6 ^ a>o(T e + T env ). 

Substituting the density Eq. 0, into the heat 

flux Q, Eq. (pQ), we obtain our first main result for the 
typical heat exchange of the Ohmic environment with 
the tunnel junction between two normal leads. The full 
temperature and voltage dependence is shown in Fig. 03 

Dynamic Coulomb interaction. — Next we discuss the 
more realistic situation where the tunneling junction is 
connected to two disordered conductors (leads). Fol- 
lowing Ref. 0|, one can find the spectral probability 
function p(oj) in Eq. (|3)) corresponding to the electron- 
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(^) 2 (2^-l)^(q,o;) 
(Pi q 2 -iuj)(V } q 2 -iuj) 



(8) 



where i,j = 1,2 are the lead indices, £>i(2) are diffusion 
coefficients within respective electrodes, and J7ij(q, oj) are 
the dynamically screened Coulomb interactions within 
(across) the electrodes. The form of spectral probability 
p(oS) \p(oS) = 2p\2 + Pw + P22] depends on the structure 
of the environmental excitations spectrum and, thus, on 
the external bias. 

The system under consideration is shown in Fig. H^i): 
two contacts are separated by distance d and their thick- 
ness is a. The external bias is V and the contacts are 
kept at temperature T and the environment at temper- 
ature Tenv Two situations are possible: i) for zero bias, 
V = we have T out = T e = T e ff = T and T in = T env . ii) 
for V 7^ the effective temperature T e ff depends on V as 
shown in Eq. (j2j). 

The screened Coulomb interaction in Eq. ([8]) in 
Fourier space has the form Uj<q,u)) — {[^^(q, + 
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FIG. 4. (color online) a): Schematic presentation of the sys- 
tem: single contact junction, with contacts consisting of two 
thin plates, which are distance d apart. Their thickness a is 
much less than the extension in x and y directions, such that 
they can be treated as 2D contacts. The temperature of the 
contacts T is kept constant, while the environment tempera- 
ture Tenv can be different, which results in heat production or 
removal in the junction, b): Heating of a tunnel junction tak- 
ing into account dynamic Coulomb interactions for the zero 
bias case (V = 0) [red lines, lower x-axis] and the voltage 
dependence for T = T e f? [green lines, upper x-axis]. The solid 
cures represent the quasi-equilibrium curves and the dashed 
assuming an equilibrium distribution for N^. (To = O.l^th, 
see text also) 
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2(q, Co>)} -1 , where ?7^(q, cj) = u(q)^+v(q)a x is the bare 
Coulomb interaction and 2(q, uj) the polarization matrix 
respectively with = V{ T>i q 2 (T>i q 2 — iu)~ 1 5iy Vi is the 
electron density of states at the Fermi surface in lead i. 

Below we concentrate on quasi 2D infinite leads. For 
this geometry with a <C where L is the characteristic 
lead size in the x and y directions, the bare Coulomb 
interaction has the form 



with z^ = (1/2 — 5n)d, leading to u(q) = 2ne 2 /q and 
v(q) = 2ne 2 e- qd /q. 

In the following, we consider the case of identical leads 
with same diffusion coefficients V\ = V2 = V and den- 
sities of states, v\ = V2 = v- The dimensionless matrix 
elements Uij of the dynamically screened Coulomb inter- 
action (in units of e 2 d) are then given by 



dzi dzj 



5(zi-z( 0) )5(z 



4 0) ) 



, (9) 



Uu 



47T 



X(<1) 



q x 2 (q) ~ coth 



X{q) coth(g) 



(10) 



where q = dq and Co 
less function x(o) = 1 



uj(d 2 / V) with the dimension- 
ing) + Al lf_ dl i x . Using these 



expressions, we can write Eq. (jSJ) as 



2e 2 d„ T 
-——ooim 
V 



I 



qdq- 



(x(q) coth(^))" 



(q 2 — iuj) 2 



(ii) 

Substituting Eq. (fTT|) into Eqs. (|3]) we can calculate the 
heat flux Q in Eq. (pQ) between environment and nano- 
j unction with dynamic Coulomb interaction. The typical 
energy scale is given by the Thouless energy for the junc- 
tion of distance <i, E t h = V/d 2 which we use to rewrite 
all expressions in dimensionless units. For a typical tem- 
perature To = 0.1£^th ~ 10 K, the temperature and volt- 
age dependence is numerically calculated and shown in 
Fig. St*). Again, the non-equilibrium heat flow Q is up to 
an order of magnitude larger and the quasi-equilibrium 
approximation Qo- We remark, that in this case the func- 
tion p(uj) introduces a natural cut-off for J(t), Eq. (|3]), 
which behaves as ~ —\t\ for large t. 

Discussion. — Above we assumed that hot electrons 
interact with acoustic phonons (acoustic environment 
modes). This assumption holds if the environment tem- 
perature T env is lower than the Debye temperature O^, 
which is of the order of optical phonon energies. In this 
temperature range the electron interaction with the en- 
vironment is quasi-elastic because the change of the elec- 
tron energy, which is equal to the energy of the emitted or 
absorbed phonons is much smaller than the electron en- 
ergy. Due to the small inelasticity of the acoustic phonon 
scattering, the deviation of the electron distribution func- 
tion in the momentum space from the isotropic one is 
small even when the electrons become hot. 



We also assume that the density of hot electrons is high 
enough so that the electron-electron scattering time r e - e 
is smaller than the time of energy relaxation r env (this 
time is large because of quasi-elastic nature of interaction 
between the electrons and environment). In this case the 
electron distribution function is close to an equilibrium 
one with an electron temperature T e , which is in high 
voltage limit is higher than the environment temperature 
T env . At very high applied voltages the electron energies 
become comparable with the energies of optical phonons 
(optical environmental modes) and the approximation of 
small inelastic scattering does not hold. 

In summary, we discussed the influence of far from 
equilibrium heating effects on properties of nanojunc- 
tions. Based on a quantum-kinetic approach we calcu- 
lated the non-linear heat flux between environment and 
junction. We showed that the resulting freezing or heat- 
ing effect far from equilibrium are by orders of magni- 
tude larger than estimates based on quasi-equilibrium 
environment theory. We obtained analytical results for 
the heat flow in an idealized high-impedance environ- 
ment and demonstrated, numerically, that these results 
hold for the more general case of an environment with 
Coulomb interaction. We showed that the environment 
can be a very effective freezing agent if the effective tem- 
perature well exceeds the high frequency cut-off ftcj max . 

One can expect that our results, in particular the gi- 
ant freezing effect, will be important for the electronic 
transport in junction arrays [9[, which will be subject of 
a forthcoming work. 
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